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The low-frequency electric microfield distribution in a Coulomb plasma is calculated for various 
plasma parameters, from weak to strong Coulomb coupling and from zero to strong electron screen- 
ing. Two methods of numerical calculations are employed: the adjustable-parameter exponential 
approximation and the Monte Carlo simulation. The results are represented by analytic fitting 
formulas suitable for applications. 
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I. INTRODUCTION 

Because of the Stark effect, stochastic electric mi- 
crofields influence optical and thermodynamic properties 
of a plasma. First, they affect the profiles of spectral 
lines and effectively lower photoionization thresholds of 
atoms and ions immersed in a plasma A compar- 

ison of experimental and theoretical widths and shapes 
of the Stark-broadened spectral lines is widely used for 
plasma diagnostics (e.g., Refs. [|[ §). Second, in some 
theoretical models of the plasma equation of state (e.g., 
Refs. H), the microfield distribution is used in order 
to calculate occupation numbers of the bound species (al- 
though such a calculation is not free from principal dif- 
ficulties, as discussed in Ref. It was shown recently 
1^ that a more accurate description of the microfields en- 
tails a considerable improvement of the equation-of-state 
model. 

In many cases, the microfield perturbation can be 
treated as quasistationary. Then the problem is re- 
duced to determination of the probability distribution of 
the low-frequency component of perturbing electric fields 
(e.g., Ref. [gl), associated with a stochastic distribution 
of perturbing ions, whereas the electrons can be assumed 
to adjust instantaneously to a configuration of the ions. 
The low-frequency microfields are appropriate to use in 
the equation of state models |^ and in calculation of 
spectroscopic line profiles for those radiative transitions 
whose frequency does not exceed the typical frequency of 
microfields produced by thermal fluctuations of the elec- 
tron density. For example, Stehle and Jacquemot used 
the model microfield method to analyze the line shapes 
and line dissolution in hydrogen plasma spectra. 
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Holtsmark [g^ has derived the microfield distribution 
function assuming that the ions are not correlated and 
the electron screening is negligible. This assumption is 
justified for very hot or rarefied plasmas, for which the 
Coulomb coupling parameter 
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is close to zero. Here, Ze is the ion charge, T is the tem- 
perature, ks is the Boltzmann constant, a = (^n^)"^/'^ 
is the ion sphere radius, is the ion number density, and 
77.20 is the electron number density (jie = Zui) in units 
of 10^° cm^"^. As we demonstrate below, the Holtsmark 
approximation is inaccurate already at T ~ 0.1. In mod- 
ern plasma experiments, T may approach unity, whereas 
in stellar matter it can be much larger. In these cases, 
correlations of plasma particles should not be neglected. 

Various approximations were developed in the past in 
order to take the ion correlations into account. If T 5 , 1, 
one may use the methods of Baranger and Mozer |l^ or 
Hooper |l|l based on a cluster expansion in powers 
of density. The electron screening is usually described by 
a Debye-like (Yukawa) effective potential, introduced in 
the context of microfield distributions by Hoffman and 
Theimer [|l^. In the limit of extremely strong coupling, 
r 3> 10, and without screening, the harmonic oscillator 
model by Mayer [ p^ is applicable, in which every ion is 
assumed to oscillate independently of the others around 
its equilibrium position at the ion-sphere center. 

The first theory capable to provide reliable numerical 
results for strongly coupled plasmas with electron screen- 
ing proved to be the adjustable-parameter exponential 
approximation (APEX) , based on a special parametriza- 
tion of the electric microfield E produced on a selected 
test particle (neutral or charged "radiator" of charge Zr) 
which undergoes the influence of charged plasma parti- 
cles ( "perturbers" of species a and of charge ZgY This 
method has been developed for Coulomb systems and 
adapted for screened Coulomb systems and ion mixtures 
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It involves non-interacting quasiparticle repre- 
sentation of the electron-screened ions, designed to yield 
the correct second moment of the microfield distribution 



(E-E) 



'Zif 
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M, (2) 



where gcr{r) and Cg. denote the radial distribution fimc- 
tion (RDF) and the relative abundance of species cr, re- 
spectively, and where kg is an effective electron screen- 
ing wave-number. After introducing the effective single- 
particle field in the form 



(3) 



the adjustable parameters {a^} are chosen to satisfy the 
condition 



4^ j e^PAPEx{e)de= {E- E). 



(4) 



The expression on the left-hand side of this equation con- 
tains the parameters {a^} to be determined, whereas 
the right-hand side can be evaluated using Eq. if 
the RDF is known. The RDF thus provides a scheme 
for evaluating the APEX microfield distribution and is 
a central ingredient whose accuracy determines the one 
of the APEX microfield results. In our implementation 
of the APEX technique, we have used the hypernetted- 
chain RDF calculations l^] . 

On the other hand, with the advent of powerful com- 
puters it is now possible to calculate the microfield dis- 
tribution from Monte Carlo (MC) or molecular-dynamics 
simulations of plasmas with the minimum of simplify- 
ing assumptions (e.g., Refs. |^ | 
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|). More- 
over, the latter methods allow one to study the effects of 
microfield nonuniformity ||2^, and to simulate high- 
frequency microfield distributions in electron-ion plasmas 
(e.g., Ref. (2^). The MC technique is based on a numeri- 
cal simulation of space configurations of a system of par- 
ticles, whereas the molecular-dynamics technique traces 
the time evolution of the system. For the low-frequency 
microfield, dynamical effects are unimportant, and the 
two methods yield the same results, as demonstrated, 
e.g., in Ref. |2^. Therefore it is sufficient to use the MC 
method in this case. 

With these powerful tools, the microfield distribution 
can be calculated now for any practically important com- 
bination of plasma parameters. However, plasma spec- 
troscopy and equation-of-state models require knowledge 
of this distribution at many different points or even in 
continuous areas of the plasma parameter space. In this 
case, either extensive numerical tables or approximate 
analytic expressions are necessary. 

We present results of calculations of the low-frequency 
microfield distribution function at a neutral and charged 
plasma point for various values of F ranging from to 



100 and for various values of an effective electron screen- 
ing length. We consider plasmas composed of a single 
species of ions; in particular, in the case of a charged 
test particle, its charge is assumed to be equal to that 
of perturbers. The calculations are performed mainly 
by the MC method; for comparison we have done also 
APEX calculations. We also present analytic formulas 
which reproduce the calculated electric microfield proba- 
bility distributions with an accuracy comparable to small 
differences between the MC and APEX results. 

In the next section, we describe basic assumptions used 
in our calculations and write down some asymptotic re- 
sults. In Sec. Ill, we present results of numerical calcu- 



lations and analytic approximations for microfield distri- 
butions produced at a neutral or charged point by ions 
interacting via unscreened or screened Coulomb poten- 
tials. The results are summarized in Sec. IV. 



II. METHOD 

A. Basic assumptions 

We consider a nonrelativistic, isotropic, overall neutral 
plasma at the thermodynamic equilibrium. The ions are 
assumed to be classical and pointlike. The electric field 
created at a point r by an ion placed at equals E{r.i — 
r) = —{Ze)~^yV{\ri—r\), where V{r) is an effective pair 
potential. This potential is taken in the Debye-Hiickel 
(or Yukawa) form: 



V{r) = {Zef 



(5) 



In the linear approximation, a test charge q embedded 
(at r — Q) creates perturbation of electron number den- 
sity he{r) = {due/ d^)ji{r), where /i(r) = — e0(r) is 
the perturbation of the electron chemical potential /x, 
and 4>{r) is the excess electrostatic potential determined 
by the Poisson equation V^0(r) = — 47r[(j(5'^(r) — eri(r)]. 
Thus in the linear (first-order perturbation) approxima- 
tion (V^ + 4TTe'^dne/d^)(j){r) — —47rqS^{r), which leads 
to the well-known (e.g., Ref. |^^) expression for the ef- 
fective screening wave number k^: 

kl ^ 4^e2 ^ = 4^ (2m,)3/2 [ksTf'^ I-,/,{x), (6) 

where Iu{x) = /o°" ^'^ dx/{e^~^ -f 1) is the Fermi integral, 
and X = f^/ksT is determined from the equation 



Ii/2{x)^7^^h^m, kB T) 



-3/2, 



=/%/2. 



The solution of Eq. (|^) and the right-hand side of Ec 
are given by accurate Pade approximations in Ref. 

In the limits of weak or strong electron degeneracy, kg 
tends to the inverse Debye length for the electrons or to 
the Thomas-Fermi wave number, respectively. 

We adopt the conventional assumption that the poten- 
tials F(|ri — r I) are additive, which is strictly valid in the 
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limit kg 0. Then the electric field E{r) is also the sum 
of elementary electric fields E{ri — r). 

It is convenient to introduce the dimensionless field /3 
and the screening parameter s: 



f3 = (a^/Ze) E, 



a kg 



(8) 



Also, /3 = (aVZe) E. 

In the canonical thermodynamic ensemble of {N + 1) 
particles, the probability density of the modulus of the 
field, (3 = |/3(ro)|, can be written as 



P(/3) = 



47r/32 
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N 



/3(ro) - /3(n - ro) 



N- 



(9) 



Here, Zn+i = / . . . / e-^(''°^''i-'''")/'=«^dri . . . drAr is 
the canonical partition function and W is the poten- 
tial energy of the configuration: W^(ro, ri, . . . , rjv) = 
1 X^i^j where Vb is the potential energy 

of the background of electrons. Our goal is to calculate 
the function P(/3) in the thermodynamic limit N ~* oo. 
In the next paragraph we shall discuss how to perform 
this calculation. 



B. Monte Carlo technique 

In the numerical MC calculations, the coordinates of 
{N + 1) particles (one test particle and N perturbers) 
are chosen in a cubic box of side length L such that 
{N + 1)/L^ — Ui. In order to include the effect of distant 
particles, the box is replicated by its "images" filling the 
space with the step L. The sum of the interaction po- 
tentials with all ion images is calculated using the Ewald 
technique and, in order to ensure isotropy, the aver- 
aging over the complete solid angle is applied. As shown 
in Ref. [^, a result of this procedure is equivalent to 
the replacement of the potential (||) in the cube by an 
effective potential Vl, 



C[,i sinh(fcsr) 
L k^r 



An 



sinh(fcsr) 



(10) 

where C^j is a numerical constant which tends to the 
Madelung constant (e.g., Ref. in the limit kg 0. 
The total potential energy W is obtained by the summa- 
tion of Vl over all pairs of particles. 

During the MC run, an ion within the box and its 
displacement are chosen randomly. If the displacement 
brings the ion outside the box, the ion is replaced by 
its image. In strongly correlated plasmas, the displace- 
ment is limited by a maximum distance smaller than L, 
in order to avoid calculation of the energy for highly im- 
probable configurations. First 10^ ion configurations are 
discarded in order to erase traces of the starting con- 
figuration. The state of equilibrium is searched by the 



Metropolis algorithm: the energy difference AW is cal- 
culated between the consecutive configurations, and the 
new configuration is accepted definitely if this difference 
is negative and accepted with probability q-^^/i^bT 
AW is positive. The latter condition allows the sys- 
tem to escape from trapping in a local energy minimum. 
When the system approaches equilibrium, W ceases to 
change appreciably. Then all equilibrium quantities de- 
pending only on particle positions (electric microfield, 
pair correlation function, etc.) can be calculated. In all 
our simulations, we try to get the maximum precision 
by minimizing the statistical errors. So we considered 
large boxes of particles (between N = 600 and N — 800) 
and, for each state (F, s), (1~6) x 10'' configurations were 
generated after equilibrium. This number of samples is 
large enough for the precision needed for all microfield 
distribution results discussed in this paper. 



C. Asymptotic and approximate expressions 

The described MC sampling procedure does not di- 
rectly provide the probability density for extremely weak 
or strong fields, which are given by rare configurations. It 
is therefore useful to know asymptotic behavior of P{(3) 
in the limits P and (3 oo. 

In the first case, the methods of Baranger and Mozer 
jnl and APEX |, |l| show a parabolic dependence of 
P{(3) near the origin. This behavior is also visible on MC 
results. Thus we assume that P{(3)/P'^ is constant near 
the origin. 

In the case of very strong fields, exact analytic results 
are available for the unscreened Coulomb potential only. 
The Holtsmark distribution, valid at F — > for any /3, 
reads Q 

P^(/3) = ^ / X eM~x^/^)sin{l3x)dx. (11) 



At /3 ^ oo, this distribution has the asymptote P{P) ~ 
1.496 /3~^/^, which is close to the asymptote of the 
nearest-neighbor (NN) field distribution 



Pnn(/3) = 1.5/3-5/2 exp(-/3-3/2). 



(12) 



In the opposite limit of extremely strong correlations 
(F — ^ oo), the Mayer model |l^ yields (for the charged 
test particle) 



Pm(/3) = v^r3/2/32 exp(-F/3V2). 



(13) 



Figure ^ illustrates the differences between various 
asymptotic theories (NN, Mayer model) and numerical 
results at finite F. Compared to the Holtsmark distri- 
bution (^l]), the most probable field values are shifted 
considerably to lower /3, the shift being much larger at 
a charged point. The NN approximation (|l^) correctly 
describes the case of large /3 for the neutral point but 
fails for the charged point or at small f3. The Mayer dis- 
tribution (O) fails to describe the high-field tail of P{0) 
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FIG. 1: Microfield distributions produced by ions interacting 
via the unscreened Coulomb potential in a plasma at F = 7 
at a neutral and charged point in the plasma. MC (solid 
lines) and APEX (short-dashed curve) numerical results are 
compared with analytic approximations by Holtsmark (long 
dashes), Mayer and nearest neighbor (dot-dashed lines). 



but provides the most probable field which, at this F, is 
offset by tens percent only. At contrast, the APEX and 
MC results are in close agreement for /3 < 2. 

The Mayer model fails in the strong-field limit, be- 
cause in this case one should consider a test ion which 
lies at a very short distance r (x /3~^/^ from the near- 
est perturbing ion. Then the geometrical and Boltz- 



mann factors give P(/3) ~ (3' 



-5/2p-r/3^ 



at (3 



In 



Appendix ^ we present der ivat ion of a more accurate 
asymptotic expression, Eq. (A4), which was previously 
given in Refs. Q for the case of the Coulomb po- 
tential without screening. A generalization of Eq. (A4) 



provides an accurate functional form of the asymptotic 
behavior of P{P) at large ^ ||, ||: 



FIG. 2: Comparison, for F = 3 and s = 0.75, of probabil- 
ity density distributions P(/3) at a charged point, calculated 
by difterent methods: Monte Carlo (solid line), APEX (dot- 
ted line), and the asymptotic expression (0) (dashed line). 
Coefficients of the asymptotic expression have been fitted on 
MC points above the cut-off value f3 — 5: F = 3.06 and 
K = 31.011. The asymptotic expression is a good representa- 
tion of the high-field tail of P{f3) since it avoids the oscillations 
shown by the two other methods (MC: numerical statistical 
noise, APEX: Fourier transform oscillations). 



III. RESULTS 

A. Coulomb potential 

In this section, we present results of MC and 
APEX calculations and analytic approximations for 
the microfield distribution brought about by the one- 
component plasma ions interacting via the Coulomb po- 
tential without screening. In this case, the probability 
density P depends on /3 and F. 



P(/3) - K /3-^/2 exp(-f - /3-3/2), (14) 

where K and F are adjustable parameters. In practice, F 
is a free parameter, whereas K can be determined from 
the normalization constraint. An example is given in 
Fig. ||. In this example, the value of F is close to the 
exact F, but this is not the general case. However, as 
shown in Refs. |2j, for any simulation of practical 
interest, it is always possible to find appropriate F and 
K to fit the high-field tail of the microfield distribution. 



1. Neutral point 

First we consider the distribution of electric microfields 
applied to a neutral test particle embedded in a plasma. 
If the plasma is weakly coupled, P{P) is given by the 
Holtsmark formula (pi]). In theoretical models (e.g., 
^) one is often interested in the cumulative probability 
distribution defined as 

/-/? 

Q(/3)= / Pi/3') dp'. (15) 
Jo 
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a, 




TABLE I: Parameters of Eq. (TO). 



FIG. 3: Microfield distributions produced at a neutral point 
by ions interacting via the Coulomb potential, for F = 0.0533, 
1, 10, and 100. Numerical results (MC, solid lines; APEX, 
short-dashed lines) are compared with analytic approximation 
(0) (dotted curves). The long-dashed curve reproduces the 
Holtsmark distribution (pT|). 



For the Holtsmark distribution, QniP) ~ 1 - 0.997/3"^/^ 
at /3 ^ (X). For arbitrary /?, accurate rational- function 
approximations to PjAB) and QniP) have been con- 
structed by Hummer [g5| . 

With increasing F, the field distribution becomes nar- 
rower, as shown in Fig. ^ The decrease of the most 
probable value pm of the dimensionless field /3, which 
corresponds to the maximum of P(/3), can be described 
by a simple approximate formula 



l-'n 



1.608 -I- 0.24 %/r 

i-H0.77Vr 



(16) 



The asymptotic behavior of P{(3) remains power-law, as 
in the case without Coulomb coupling. This facilitates 
construction of self-consistent rational approximations to 
(5(/3) and P(/3). We have calculated P{P) for various F 
from to 10 by the APEX method and for F up to 100 
by the MC technique. Our fit to Q{P) reads 



QiP) 



go 13^ - 1.33 P^/^ + 
+ q2P^ + (73 ^3 - + pe 



(17) 



where (?„ — an{l + Pn •\/r)~^", and the parameters a„, 
Pm and 7„ are given in Table |. P{P) is obtained from 
Eq. ( |l7| ) by elementary differentiation. At F = 0, this 
differentiation reproduces Ph{P) at any p with a maxi- 
mum fractional error of 0.24%. At finite F, the difference 
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FIG. 4: Microfield distributions produced at a charged test 
particle by ions interacting via the Coulomb potential at var- 
ious r from 0.1 to 100. MC results (solid lines) are compared 
with analytic approximation (|o|) (dotted curves). The Holts- 
mark distribution (F — 0) is also plotted (dashed line). 



between the fit and the MC data increases up to several 
percent, remaining however not larger than the differ- 
ence between the MC and APEX results, as one can see 
in Fig. |. 



2. Charged point 

For the charged test particle, the asymptotic behavior 
of P{P) at large P is qualitatively different in the cases 
of zero and non-zero F. When F 7^ 0, the power-law 
decrease of P{P) is replaced by an exponential. At mod- 
erate P and large F, P{P) is approximately described by 
the Mayer distribution (O). The cumulative function of 
this distribution is 



Qm(/3, r) = erf (/3 ^172) - ^ 



Pe 



(18) 
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FIG. 5: The same as in Fig. W but on the logarithmic scale. 



which is easily calculated using, e.g., the highly accu- 
rate rational approximation to erfc(a;) in Ref. 
The most probable value of the microfield provided by 



this distribution is (3^ = a/Z/T. With increasing F, /3m 
for a charged point decreases faster than (3^^ given by 
Eq. The MC results for f3m are approximately de- 

scribed by the following modification of f3j^: 



f3^J^ = Teff^ 0.774 + ri/4 + r. (19) 

The microfield probability density at various values of 
r is shown in Fig. ^. The exponential decrease of P{/3) 
at large F is more easily seen in the logarithmic scale 
(Fig. 1^). According to Eq. (|l^), the exponent at large 
f3 is pro portional to and not to as in Eqs. ( p^ ) 
and ([l8|). Nevertheless, it is still possible to construct a 
self-consistent analytic approximation to Q(/3) and P{(3) 
analogous to Eq. (|l^) . Since the asymptotes are qualita- 
tively different in different coupling regimes, this approx- 
imation is more complicated: 



Qo(/?) + 0.873 VrQA/(/3,Feff) 
1-K 0.873 \/r 



(20) 



where 



m = 



1/3^ exp(-F'/3i/2) +/36 



2.257rg (1 + F0-6)-2-^5 + 15.3 f3^ + 1.238 q/S^ + f3^/^ exp(-F'/3i/2) + 

I 



(21) 



g = 9.19 -f 2.178 F^-^'', and F' F/(l + 0.19 F°-^^^). applications; the results are presented below. 



At F = 0, dQa/dp reproduces Ph{P) within 1%. The 
accuracy deteriorates at F sa 0.1 but improves again at 
higher F: in Figs. ^and|^, differences between the fit and 
MC results are barely visible at F = 0.5, 3, 50, and 100. 



B. Effective screened potential 

In this section, we consider microfield distribution in 
an electron-screened Coulomb plasma. Assuming that 
the ions interact via the effective potential (^), we have 
performed MC simulations of P{(3) for various values of 
the screening parameter from s = 0.05 to s = 3.0 and the 
Coulomb coupling parameter from F = 0.1 to F = 100. 
In accord with an intuitive expectation, the MC simula- 
tions show that the typical fields applied to a test par- 
ticle are reduced when the electron screening is taken 
into account. In this case, the probability density P and 
its cumulative function Q depend on three dimensionless 
arguments: /3, F, and s. Naturally, analytic approxi- 
mations in this three-dimensional space become complex 
and less accurate than at fixed s = 0; nevertheless we 
have attempted to construct unified formulas for evalua- 
tion of P{I3) with an accuracy which is sufficient for most 



1. Neutral point 

The most probable field strength applied to a neutral 
test particle, evaluated by the MC method, can be pa- 
rameterized as 



/3r(r,s)«/3i°)(s) 



1 + o.i5\/r 



H-0.77(l + s)e-i-s« Vr' 



where 



/3m' (s) ~ [0.622 -f- 0.25 se"]" 



(22) 



(23) 



Equation ( |22D extends Eq. (g^ to the case where s ^ 
0. It is valid for the whole considered range of plasma 
parameters, F < 100 and s < 3. 

The reduction of the typical microfield strength applied 
to a charged test particle, when the screening is taken 
into account, is illustrated in Fig. ^, where the dashed 
lines correspond to the case of a neutral test particle, 
which we consider in this section. At a constant F, (3m is a 
monotonically decreasing function of s. The dependence 
of Pm on F is less obvious. At small values of s, f3m 
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FIG. 6: The "most probable" field strength [normalized ac- 
cording to Eq. ^] as a function of the Coulomb coupling pa- 
rameter r at four values of the screening parameter s. Dashed 
lines: neutral point; solid lines: charged point. 
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FIG. 8: Microfield distributions produced at a neutral point 
for the coupling parameter F — 0.5, for 7 values of the screen- 
ing parameter s marked near the curves. MC results (solid 
lines) are compared with the analytic approximation (dotted 
lines) . 
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FIG. 7: Microfield distributions produced at a neutral point 
in a plasma by ions interacting via the screened Coulomb 
potential for two values of the screening parameter s (dashed 
lines: s = 0.05; solid lines: s = 1.5) and two values of the 
Coulomb coupling parameter F. The analytic approximation 
(Appendix S) is shown by dotted lines. 



decreases monotonically with increasing F. However, the 
opposite is observed when s is large: in this case, /3m 
increases with growing F. This impHes that the most 
probable field strength depends on s stronger at small F 
and weaker at large F. 

The modification of the probability density profile with 
variation of the plasma parameters F and s is shown in 
Figs. |7|-|ic|. In the most important range of the coupling 
parameter, F < 10, the dependence of P{(3) on F becomes 
slow, as the screening becomes sufficiently strong. For 
example, in Fig. |^ we observe a significant modification of 
P{P) at s = 0.05, when F increases from 0.1 to 5, whereas 
P{P) is only slightly modified with the same increase of 
F, if s = 1.5. 

On the other hand, the profile of P{(3) strongly de- 
pends on the value of the screening parameter s (Fig. ||) , 
especially if s > 1. In this case, the most probable field 
is reduced drastically; simultaneously, the distribution 
acquires a long "tail," which shows that the values of 
(3 ^ (3m occur more often than in the s = case. Fig- 
ure ^, which presents the same dependences on the loga- 
rithmic scale, clearly reveals the two limiting power laws, 
P (X (3^ and (3~^^'^ at small an d la rge (3 values, respec- 
tively, in agreement with Sec. |II C| . We see that these 
limits are approached considerably more slowly in the 
case of strong screening. 

A comparison of Figs || and |l^ shows that P is less 
sensitive to s, when F is large, in accord with the afore- 
mentioned property of (3m- Nevertheless, even at F 3> 1, 
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-0.5 
logio /5 




0.5 



FIG. 9: The same as in Fig. M on the logarithmic scale. 




-0.5 
loglo ^ 

FIG. 10: The same as in Fig. | but for F = 50. 



the sensitivity of P with respect to s remains essential. 

It is possible to construct a fitting formula to the prob- 
ability function Q{(3) analogous to Eq. (p^), taking into 
account the screening. Such parametric approximation 
is given in Appendix ^ and is shown in Figs. [7]-pX| by 
dotted lines. Although less accurate than the fits pre- 



FIG. 11: Microfield distributions produced at a charged test 
particle by ions interacting via the unscreened (dashed lines) 
and screened Coulomb potential (solid lines; the screening 
parameter s = 0.7), for F — O.L 1, and 5. The analytic 
approximation given in Appendix O is also shown for all cases 
(dotted curves). 



the numerical results obtained by the MC simulation. 



2. Charged point 

Consider now the microfield distribution created by 
plasma ions at a point where one of these ions is placed, 
assuming that the Coulomb interaction is screened ac- 
cording to Eq. (|^). As well as in the case of a neutral 
point, the screening lowers typical microfield values. De- 
pendence of the most probable field (3m on s and F can 
be approximated by a simple expression 



/3™(r,s) = /3}„"U.) 



(0), 



where /3^'(s) is given by Eq. (^. At s 0, Eq. (|2j) 
reproduces Eq. (n9|) . The dependence of (3m on F at var- 
ious values of s is plotted in Fig. |^ by solid lines. Unlike 
the case of a neutral point considered in the preceding 
section, (3^ decreases monotonically with increasing F 
at any given value of s. 

The modification of the probability density P{P) with 
increasing s or F is illustrated by Figs. pl]-p^. The depen- 
dence on the coupling parameter F remains qualitatively 
the same as without screening: with the increase of F, 
the typical field strengths become lower, and the distri- 
bution P{(3) becomes narrower. The increase of s also 



ri/4 , 



0.774 -h 0.54 se« 



-1/2 



(24) 



sented in Sect. Ill A, this approximation reproduces well 
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FIG. 12: Microfield distributions produced at a charged test 
particle for the coupling parameter F = 10, for the values 
of the screening parameter s, marked near the curves. MC 
results (solid lines) are coinpared with the analytic approxi- 
mation given in Appendix O (dotted lines). 



FIG. 13: Comparison of distributions P{f3) at two values 
of r (solid lines: F = 0.3; dashed lines: F = 10) and seven 
values of s (s = 0.05, 0.50, 1.04, 1.5, 2.0, 2.5, and 3.0, from 
right to left) on the logarithmic scale. Dotted lines show the 
approximation (Appendix |c|) . 



shifts the peak of P to smaller /3. At s > 1, however, the 
latter shift is accompanied by a striking modification of 
the shape of the function P: a fast growth at small (3 is 
followed by a slow, gradual decrease at /3 > /?,„. 

As in the case of the Coulomb potential, the limiting 
behavior of P(/3) at /3 ^ changes from power law 
at r = to the exponential decrease oc p~^/'^e~^^ at 
/3 — > oo. This limiting law is reached very slowly, if s is 
large, as clearly seen in the logarithmic scale (Fig. 13). 
When f3 is moderately large and F ^ 1, the decrease is 
approximately Gaussian, cx e~^^ 

The rich variety of the shapes of P{f3) (depending on 
s and F) complicates significantly the construction of a 
fitting formula. In this unified fit to Q and P (like 

those presented in the previous sections) would become 
too cumbersome. Therefore we have chosen to construct 
an analytic approximation to the function P only. When- 
ever necessary, Q can be found by numerical integration 
[Eq. (HI)]. 

Our approximate formula is presented in Appendix |^. 
Its quality is revealed by Figs. |l^-|lj, where the fit is 
compared with results of the MC simulations. The typi- 
cal accuracy of several percent at s < 1.5 is expected to 
be sufficient for most applications. At s > 1.5 the accu- 
racy deteriorates, and at s > 2 the asymptotic behavior 
at large /? is not reproduced, as one can see from the log- 
arithmic plots (Fig. |l^. At such strong screening, the 
fit still may be used for evaluation of P{f3) not too far 



from Pm- Indeed, on the linear scale (Figs. [Ij) the 
difference between the fit and the MC results appears to 
be small even at s = 3. 

Figure |l^ allows one to compare the screening effects 
in the two cases of a neutral and charged test particle. 
Since F is greater than unity, the difference between P{(3) 
functions in the two cases is large at s = 0.05 (i.e., for a 
nearly Coulomb potential), in agreement with Sec. Ill A. 
With the increase of s, however, the difference becomes 
smaller near the peak of P{0)- The positions of the peaks 
for a neutral and charged points almost coincide at s > 
1, and the difference in their heights is caused by the 
fact that, in the case of a neutral point, P{(3) decreases 
much more slowly at large /3, and therefore the region of 
P 3> Pm gives a larger contribution to the normalization 
integral, than it does in the case of a charged point. 



IV. SUMMARY 

We have calculated microfield distributions at neutral 
and charged test particles in a one-component plasma 
of ions, interacting via Coulomb potential, in various 
regimes from weak to strong coupling. The MC and 
APEX methods of calculation yield similar distributions, 
in agreement with previously known results Self- 
consistent elementary-function approximations for the 
field probability density P(/3) and its cumulative distri- 
bution Q(/3) are constructed in the two cases of a neutral 
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charged 
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fit 




dimensionless screening parameter s [Eq. (Q)] varies from 
to 3. The whole set of numerical results for P(/3) at 
various values of the coupling and screening parameters 
is approximated by analytic expressions. 

The obtained results can be used in theoretical mod- 
els of optical spectra and equations of state of Coulomb 
plasmas. 
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APPENDIX A: HIGH-FIELD ASYMPTOTIC 
EXPRESSION FOR THE MICROFIELD 
DISTRIBUTION 



FIG. 14: Comparison of the probability density distributions 
P(/3) at a neutral (dashed lines) and charged (solid lines) 
point at r = 1.5 and for four values of s. Dotted lines show 
the approximations in Appendixes ^ and |c|. 



and charged point, for a Coulomb coupling parameter F 
varying from to 10^. 

Furthermore, MC calculations of the microfield distri- 
bution have been performed for the screened Coulomb 
interaction, using the model of ions interacting via the 
Debye-like (Yukawa) effective potential, with an effective 
screening length as a second independent parameter. The 



In this Appendix, we outline the derivation of the 
strong- field asymptotic limit of Eq. (^). For brevity, we 
adopt the convention that all lengths and radius vectors 
are measured in units of the ion-sphere radius a. 

We take advantage of the fact that, for Coulomb inter- 
actions, high-field contributions are produced by nearest 
particles. This well-known result has been investigated 
in Rcf. ||l^ . We assume (i) that the microfield at rp = 
is dominated by the contribution of the nearest neighbor 
ion located at ri, and (ii) that only the potential of this 
ion contributes to the potential energy. Then Eq. (||) can 
be approximated by 



J 



P(/3) ~Pa«(/3) = 



^ fo dri f ^ uro . 

2 JO ^ 'fr2>ri ^ 



dro 



fn dri f ^ dr2 ■ ■ ■ f 



Ir2>ri 

r 



dvN e-r'"('-i) 



(Al) 



Here, D, is the total volume of the system, and v{r) = 
V{r) a/{Ze)^ is the reduced potential, so that Tv{r) = 
V[r)/kBT. In our units, n = ^ttN. 

Taking into account that j^dr ~ il, 

n{l-rl/N), and limjv^oo(l-r?/A^)^ = 



i'2>ri 



dro = 



we obtain 



In 



5{(3-Vv{r))dr 



Jo 



-rv{r) 



(A2) 



The coordinate transformation u — Vti(|r|) 
(3{r)r/r with the Jacobian J = {r"^ / P'^)\dr / d(3\ yields 



1 

K 



dr 
dp 



-rv{r) 



(A3) 



with /3 = /3(r) = \dv{r)/dr \ and K ^ ^2g-r3-ri,(r) 

For a Coulomb potential, v{r) — 1/r and /3 = l/r^, so 
that 

/3-5/2 exp(-F/3i/2 „ ^-3/2) 



2/„°°/3'-^/2 exp(-r/3^ 



,,1/2 



For a Yukawa potential v{r) = e 



P{P) 



1 



exp[ 



Fe- 



- l3'-^''^)d(3' 
"/r, we have 
'7r] 



K {2 + sr + 2/sr) {s/r'^) e""'' ' 
where r should be determined from the equation 

/3=(l + sr)^, 



(A4) 



(A5a) 



(A5b) 
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and 



K = 



' cxp[— - Te *'"/r] dr. 



(A5c) 



Equa tion ( [A5| ) is considerably more complicated than 
Eg. (A4); it can be compared to Hoope r's f ormulation 
||13| . On the other hand, the simpler Eq. (A4) or its gen- 
eralization, Eq. (]r^), can be sufhciently accurate for most 
applications. This has been verified by comparing with 
MC results pi, Q and is illustrated in Fig. |. 



APPENDIX B: APPROXIMATION TO THE 
PROBABILITY FUNCTION OF MICROFIELDS 
AT A NEUTRAL POINT 

In this Appendix we present a fitting formula to the 
probability function Q{/3), Eq. (p^), for a neutral point 
in a plasma with Coulomb coupling and screening. 

At every pair of T and s values, we derive a Fade ap- 
proximation to the microfield probability function. 



Q(/3) 



ai + 02/3 + 03/32 + a4/33 - P^/^ + (3^ 

This expression ensures that, when its derivative is taken, 
the Hmits P{(3) cx /J^ at ^ ^ and P(/3) « l.S/?-^/^ 
(3 ^ 00 are reproduced. At arbitrary /?, an agreement 
with MC results is provided by an appropriate choice of 
the fitting parameters 09-04. The latter parameters, in 
turn, can be approximated as functions of F and s: 

+ (59 + 8.1,s2)g, (B2) 



Oo 



oi 



a2 



as 



04 



97 s' 



1 + 3.1 X 10- 



1.16 



a{s) 



0.188 s6 
0.068 



1035' 



4s) 



l-f 0.33 s 
0.038 



1-1- 0.030 



95 s 



1 + 6 X 10-3 
1.894 -Ks 



1.2 s^g. 



2 -t- s 



(B3) 

(B4) 
(B5) 
(B6) 



This approximation has been checked for the whole 
range of the plasma parameters for which the MC simula- 
tions were performed, i.e., at < F < 100 and < s < 3. 
In the case of purely Coulomb potential (s = 0), however, 
Eq. ( p^ ) should be used as more accurate. 



APPENDIX C: APPROXIMATION TO THE 
PROBABILITY DENSITY OF MICROFIELDS AT 
A PLASMA ION 

In this Appendix we present an analytic approximation 
to the probability density P(/3) of electric microfield at 
an ion in an electron-screened Coulomb plasma. 



At every F and s, we write 



'AT 



Be 



-bp'< 



-VP 



1/2 



l + c/39/2 



, (CI) 



where Sn is the normalization constant. For the latter 
constant, we have 



S^=A + B -f T~^F{c/T% (C2) 



aa 



3/q 



^53/7 



where F(3/a) and r(3/7) are the Gamma-function values 
which are easily calculated (e.g., Rcf. p^), and 



F{y) 



00 3.2 g-v^ 
1+2/X9/2 



dx. 



(C3) 



For the latter integral, we have constructed an approxi- 
mation. 



(Bl) F{y) « 1 



47r 



,1/9 



1 

240 



0.849 y^/^-f 3.2 y^/^+ 2.43 y^'-^ y 



7/9 



• (C4) 



Expression (|CJ) fits the integral (|C^) within 0.8%. This 
is sufficient for evaluation of Q(/3) and P{P) in most ap- 
plications. The accuracy of Eq. ( |C4| ) may be insufficient, 
however, if the values of [1 — Q(/3)] at /3 ^ /?„ are of 
interest. In this case, the normalization c ons tant can be 
evaluated numerically according to Eqs. (C2) and (|C3|). 

Equation ( |Cl|) ensures that P{[3) cx /J^ at /3 ^ 0. 
Moreover, it also ensures the correct limiting behavior 
at /3 — > 00, Eq. (|lj), with f = F, provided that a and 
7 are both greater than 0.5. With the choice of param- 
eters presented below, this is the case at < s < 2.3 
(any F), which covers the whole range of values typically 
encountered in stellar and laboratory dense plasmas. 

We have parameterized A, a, a, B, b, and 7 in Eq. (CI) 
as functions of F, having the same form at any s: 



i + AiVr 

A = A-i : ^ : r. 

1 + A2F2 + ^3F4 

o = oo + F/2, 
ai + 2a2 

^ ^ 1 + B2 F2 + S3 F4 ' 
b = 5o + F/4, 

_ 71 + 1.572 \/r 
^ " 1 + 72 \/r 



(C5) 
(C6) 
(C7) 

(C8) 
(C9) 

(CIO) 



Then the parameters of these expressions (A„, _B„, oq, 
60, ai, a2, 71, and 72) and the parameter c of Eq. (CI) 
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have been approximated by analytic functions of s: 

Ai = 0.59 + 2540 5^ + 3 5^'', (Clla) 

.^5 -I- 1 n -J- 9 .<i4-5 
A2 = - 



^3 



l + 20sO-5 
2.17 X 10-3 s^ 
= 14.8/[1 + 117s3-5], 

ao = 1.15 + 2 s^^ 



(Cllb) 

(Cllc) 
(Clld) 



(C12) 



bo 



(1 + 0.54 s^-^)/(l + 0.07s), 



71 = 0.1 + 1.1/(1 + 0.174s^-5), 

5.4 1.1 

72 = 



1 + 21si 



1 + 19s' 



0.16 ' 



(C15) 



(C16a) 
(C16b) 



ai 
a2 



B2 



0.1 + 1.1/(1 + 0.145 s^), 
5.4 1.1 

1 + 20s2 ^ 1 + 14s0-35' 



0.386 + 300 s2 + 1.1 s^•^ 
0.038 + 0.79 s°■^^ 

3.7 X 10-3 s5.5 



1 + 4 X 10-3 s9' 



(C13a) 
(C13b) 



(C14a) 
(C14b) 

(C14c) 



0.097 



l + 210s2-5 exp(-1.3si- 



(C17) 



The high powers of s in some of these equations effec- 
tively describe the strong s-dependence of the shape of 
the function P{(3) at s > 1. 



As noted in Sec. IIIB2, this approximation is valid at 



r < 100 and s < 2, but it can be also used at 2 < s < 3, 
provided that /? is not larger than « 10/3^. 
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